Density correlations and dynamical Casimir emission of Bogoliubov phonons in 
modulated atomic Bose-Einstein condensates 
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We present a theory of the density correlations that appear in an atomic Bose-Einstein condensate 
as a consequence of the dynamical Casimir emission of pairs of Bogoliubov phonons when the atom- 
atom scattering length is modulated in time. Different regimes as a function of the temporal shape 
of the modulation are identified and a simple physical picture of the phenomenon is discussed. 
Analytical expressions for the density correlation function are provided for the most significant 
limiting cases. This theory is able to explain some unexpected features recently observed in numerical 
calculations of Hawking radiation from analog black holes. 
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The dynamical Casimir effect [l| is a very general pre- 
diction of quantum field theory: whenever the boundary 
conditions and/or the dispersion law and/or the back- 
ground of a quantum field are quickly varied in time, 
pairs of quanta are generated off the vacuum state by 
parametric amplification of zero-point noise. The sim- 
plest and most celebrated example of dynamical Casimir 
effect was predicted for an optical cavity whose plane- 
parallel mirrors are made to rapidly oscillate in time 
along the cavity axis Hi, 3. Despite the significant ef- 
fort devoted to these fascinating effects, no experimental 
observation of the dynamical Casimir effect has yet ap- 
peared, the main reason being the difficulty in moving the 
mirror at a fast enough speed 0]. Alternative schemes 
to modulate the effective optical length of a cavity on 
a very fast time-scale by acting on the refractive index 
of the cavity material have been proposed and shown 
to give a sizeable intensity of dynamical Casimir emis- 
sion 

S I, 0, 1 ©, EE E3 ■ Experiments in this direction 
are in progress in several groups [H, EH ■ 

Since the original proposal by Unruh [f4j], the ad- 
vances in the field of the so-called analog models [f5[ have 
pointed out the possibility of simulating the physics of a 
quantum field on a generic curved space-time in table- 
top condensed-matter experiments: the propagation of 
elementary excitations in spatially and temporally inho- 
mogeneous systems can in fact be recast in terms of a 
relativistic wave equation on an effective curved space- 
time. In the simplest case of acoustic waves in a fluid, 
the space-time metric is fixed by the spatial and temporal 
profiles of the sound speed and of the flow velocity. Upon 
quantization of the resulting field theory, an analog dy- 
namical Casimir effect is then expected to appear when- 
ever the sound speed in a spatially homogeneous system 
is made to quickly vary in time, which in the language 
of the analogy corresponds to the expansion or contrac- 
tion of the underlying universe 0, E3 ■ On the other 
hand, in the presence of an acoustic horizon separating 
an upstream region of sub-sonic flow from a downstream 
one of super-sonic flow, the emission of analog Hawking 



radiation has been predicted [3, EH, EH • 

As proposed in 0, [ll| , a most promising way of exper- 
imentally investigating this physics involves the measure- 
ment of the correlation function of density fluctuations. 
A recent numerical experiment [22| has observed both 
the dynamical Casimir effect and the Hawking emission 
in a microscopic simulation of the dynamics of an atomic 
Bose Einstein condensate during and after the formation 
of an acoustic black hole. 

After a scries of papers [lj| E3, HI, [H| focussed on 
the deposited energy and the spectrum of the emitted 
phonons, a number of theoretical works have started in- 
vestigating the time-evolution of the correlation functions 
of atomic many-body systems after a sudden quench of 
some parameter, in particular the atom-atom interaction 
constant [2(| HH, Hy, [27|, 28]. In the meanwhile, several 
experiments have recently addressed the response of Bose 
gases to a time-modulation of the atom-atom scattering 
length: the interest of most of these experiments was 
however concentrated on the energy deposited in the sys- 
tem by the perturbation [2^] and only provided qualita- 
tive information on the created density perturbation [3(| . 

As density correlations are becoming a standard ob- 
servable in the experimental study of ultracold atomic 
gases [3l], HI, HH we expect that the measurement of the 
density correlation pattern created by suitable perturba- 
tion sequences will soon provide a new powerful tool to 
investigate the microscopic properties of strongly corre- 
lated atomic gases and in particular of their elementary 
excitations. A few pioneering studies in this direction 
for the simplest case of a sudden quench have recently 
appeared p7l. [28j. 

It is remarkable to note that a similar strategy is 
presently under way in the completely different field of 
cosmology, where one is trying to extract information on 
the primordial inflation phase of the universe from the 
spectrum and the correlation function of the temperature 
fluctuations of the cosmological microwave radiation [34| . 
In the language of the analogy, the fast expansion of the 
universe in the inflation phase corresponds in fact to a 
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fast temporal variation of the sound speed, which is in- 
deed the result of a fast variation of the atom-atom scat- 
tering length [IE US, 13 . 

The present paper presents a comprehensive theoret- 
ical investigation of the different features that appear 
in the density correlation funcion of a spatially homo- 
geneous Bose Einstein condensate as a consequence of a 
time-dependent atomic scattering length. Calculations 
based the Bogoliubov theory of dilute condensates allow 
for a quantitative understanding of the different regimes 
as a function of the temporal profile of the scattering 
length modulation and provide a clear physical picture 
of the phenomenon in terms of the dynamical Casimir 
emission of entangled pairs of phonons. This theory fully 
confirms the numerical observations that was put forward 
in the original paper [22j ■ 

The paper is organized as follows. In sec|l]we present 
the physical system under investigation and we briefly re- 
view the Bogoliubov approximation. The general theory 
of the density correlations that result from the dynami- 
cal Casimir emission is presented in Sec|TTJ The following 
sections discuss the phenomenology in the most remark- 
able cases of an adiabatic transition (Sec lIIIj) . of a sud- 
den jump (Scc lIVj) and of a slow ramp (Scc|V| of the 
atom-atom scattering length. Analytical formulas valid 
in the hydrodynamic limit are presented in Sec lVIl A 
comparison with the numerical results of [22I ] is presented 
in Sec lVIII The case of a quasi-periodic modulation of 
the scattering length is analyzed in Sec I Villi and a pos- 
sible application to the measurement of temperature dis- 
cussed. The possibility of reinforcing the density corre- 
lation signal by mapping phase fluctuations into density 
ones is quantitatively studied in Sec lIXl Conclusions are 
finally drawn in SecfXl 



by the Gross-Pitaevskii equation, 

..<9</>o h 2 V 2 , Tr . ^h 2 a{t). ± . 2± M . 

*fi-5T = --5 — fo + V(x)<j) + l±\fo\2fo. (1) 

at 2m to 

In the spatially homogeneous case V(x) = that we are 
considering here, the condensate wave function remains 
at all times constant in space and only acquires a time- 
dependent global phase, 0o(x, t) = \fn expi8(t). 

Within the Bogoliubov approximation [37| , the fluctu- 
ations around the purely condensed state are described 
by a quadratic Hamiltonian in the Ak operators describ- 
ing the non-condensed k ^ plane-wave modes 



H = E + 



(A k ,A_ k ) 



£ k + fi(t) /z(i) 
-fi(t) -E k - 



A k 
At 



(2) 



The A k operators have a simple expression in terms 
of the momentum-space atomic field operators, Ak = 
A -1 / 2 & k =o bk' Here, N is the total number of particles in 
the gas. In the following of the paper, we shall also make 
use of the real-space operators A(x) that are defined as 
the Fourier transform of the Ak operators (3?| • The (in- 
stantaneous) chemical potential is defined in terms of 
the nonlinear interaction constant g(t) — 4irh 2 a(t)/m as 
fi(t) = g(t) n. The kinetic energy of the k-mode is indi- 
cated as E k = h 2 k 2 /2m. 

Neglecting the zero-point energy, the Hamiltonian ([2|) 
can be recast for each instantaneous value of a(t) into the 
canonical form 



(3) 



I. THE PHYSICAL SYSTEM AND THE 
BOGOLIUBOV DESCRIPTION 

We consider a spatially homogeneous Bosc-Einstein 
condensate of atoms of mass to. The gas is assumed 
to be initially at rest in its thermal equilibrium state at 
a temperature T and to have a density n. Atom-atom 
interactions are modeled by a repulsive, local interac- 
tion potential of scattering length a > 0. The value of 
the scattering length is assumed to be constant in space 
but to have a non-trivial temporal dependence a(t). In 
the remote past and future t = ^00, it tends to a con- 
stant values. Experimentally, the possibility of tuning 
of the atom-atom interactions on a wide range has been 
demonstrated using magnetic and optical Feshbach reso- 
nances [3(| , as well as by modulating the lateral confine- 
ment of reduced dimensionality samples [3l| . 

At all times, the system is assumed to be well within 
the dilute regime noa 3 <C 1 where the time-evolution of 
the condensate is accurately described by the so-called 
Bogoliubov approximation |3ll . |37| . The time-evolution 
of the classical condensate wavefunction </>o is described 



where the (bosonic) Bogoliubov operators a k (a k ) respec- 
tively destroy (create) an elementary Bogoliubov excita- 
tion of momentum ftk and energy 



huk{t) = ^/E k (E k + 2fi(t)). 



(4) 



In terms of the atomic field operators, the (instanta- 
neous) Bogoliubov operators Sk have the following ex- 
pression: 



a k = u k (t) A k - Uk(t) Al 



(5) 



in terms of the (instantaneous) Bogoliubov coefficients 
«k(<) and Vk(t), 



±1/2 



(6) 



as a consequence of the time-dependence of the scattering 
length a(t), the Bogoliubov operators have an explicit 
time-dependence even in the Schrodinger picture of 
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At each time, the (instantaneous) ground state \g(t)) of 
the Bogoliubov theory is defined by 



a k \g) =0 Vk. 



(7) 



The time-dependence of \g(t)) is the key responsible for 
the dynamical Casimir emission: when the scattering 
length is modulated at a fast rate, the system is not able 
to adiabatically follow the istantaneous ground state. 
Non-adiabatic processes then result in the creation of 
correlated pairs of excitations in the system out of the 
vacuum state This point of view will be discussed in 
full detail in SccflVl 

An alternative, yet equivalent picture of the dynamical 
Casimir emission can be obtained in the limiting case 
of a weak modulation a(t) — + Sa(t) with |5a(i)| -C 
a . In analogy to the discussion of Q, the dynamical 
Casimir emission can in this case be recast in terms of 
the following Hamiltonian 



Ti = Tio 



2irh 2 n . , 2 
6a(t) > (u k +Vk) 2 



k 



x (a k + a_ k )(a k + a^ k ). (8) 



Here, Ho is the Hamiltonian of the gas for a constant 
value of the scattering length cto and the effect of the 
weak modulation Sa(t) is to simultaneously excite pairs 
of entangled Bogoliubov particles with opposite momenta 
±k. In the language of nonlinear optics, such a process 
goes under the name of parametric down-conversion [44| . 
The most remarkable case of a periodic modulation of 
a(t) will be discussed in Sec lVIIIl 



II. THE THEORETICAL FRAMEWORK 



A. The time-evolution of Bogoliubov operators 



via the time-dependence of the w k and i; k Bogoliubov co- 
efficients, and are responsible for the mixing of the a k 
and a^L k operators. These mixing terms are proportional 
to the rate at which a(t) is varied in time. Once a(t) has 
approached its late-time limiting value, one is left with 
the trivial oscillation of the a k (£) and a_ k (i) operators 
at frequencies ±w k . 

The relation between the a k (£) and al k (t) operators at 

a generic time t to their initial values a k (0) and a^ k (0) 
before the excitation sequence can be summarized as a 
pair of linear equations. 

&k(t) = [c k ,+(t)a k (o) + c k ,_(t)6i k (o)] (ii) 

= fe + (t)4(0) + Pk,-(*)6-k(0)l, (12) 



whose coefficients Ck,±(t) have to be computed by solv- 
ing the pair of differential equations ([MTU]) . 

Thanks to the thermal equilibrium hypothesis, the Bo- 
goliubov modes are assumed to be initially uncorrelated 
and thermally occupied, 



<a k (0)a_ k (0)) = 
<4(0K(0)) = r£'° 



1 



cxp{hwk/ UbT) — 1 



(13) 
(14) 



At late times t > t& n after the end of the modulation, 
the Bogoliubov mode occupations 



ra k (t) 



(o£(*)3k(f)) = 

= (|C k , + (^ n )| 2 + |C k: _(t fin )| 2 ) 4 h '°+ 

+ |Ck,-(<fin)| 2 (15) 



Thanks to the quadratic nature of the Bogoliubov 
Hamiltonian ||3J), the time-evolution of the Bogoliubov 
operators in the Heisenberg picture can be written in a 
closed form: 

&"k . / . . \ *+ 

= -iw k a k + (Ukv k -u k v k ) a'_ u . (9) 



da 1 



dt 



= iw_ k ol k + (uk «k - u k Uk) Ok- (10) 



At each time t, the zt k (£) and w k (t) functions are to be 
evaluated using the instantaneous value of the scattering 
length a(t) according to (J3J). Dots indicate the derivative 
over time t. 

The first terms on the RHS of (f9"l [TTj|) describe the triv- 
ial evolution of the a k and a^L k operators at frequencies 
±o; k under the instantaneous Hamiltonian ||3J). The other 
terms take into account the dependence (O of the a k and 
a^ k operators on the instantaneous scattering length a(t) 



remain constant in time, while the anomalous averages 



A k (t) = (a k (t) o_ k (i)> = C k ,+(t fin )C k ,_(i fin ) 



x (2n 



th,0 



1 e 



-2lUJ\^ (t — tfi n ) 



keep on oscillating at the frequency 2cjk- 



(16) 



B. The density correlation function 

In a homogeneous system the modulation of a(t) has 
no effect on the average density that remains flat at all 
times, 



i(x) = (* t (x)*(x)) = n. 



(17) 
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On the other hand, the emission of entangled phonon 
pairs is clearly visible in the density correlation function 



5 (2) (x,x') = 4 /# t (x)^ t (x')*(x')*(x: 

= 1 + ^E[ eik(x " x) (K+- k ) 2 + i) 



1 

'nV 



y^(Mk+^kj" 



a k a k ; 



x ^(aj c + a-k)(a k + aL k )^> 



1 



(18) 



which involves a combination of Bogoliubov mode occu- 
pation (p~5|) and anomalous average (fT6|) . As we shall see 
in the following of the paper, the different time depen- 
dence of the two terms is responsible for qualitatively 
different features in the density correlation function. 



C. Effect of external trapping 

Before proceeding with the analysis of the density cor- 
relation function, it is important to briefly clarify the 
consequences of the spatial inhomogeneity of the gas in 
the presence of a trapping potential. 

As the density profile of a trapped gas strongly de- 
pends on interactions, the modulation of the scattering 
length a(t) may result into a macroscopic oscillation of 
the condensate shape [38| and even in its collapse when 
the scattering is tuned to a large and attractive value. 
This latter effect has been experimentally demonstrated 
in a remarkable way in the so-called Bose-nova experi- 
ments of nn. 

In the framework of the standard Bogoliubov the- 
ory [3^ | , this physics is described by a term of the form 



SHi 



2irh 2 Sa(t) 



d J x|0o(x)|^5(x)A(x)+h.c (19) 



where </>o(x) is the condensate wavefunction, normalized 
in a way that Jd 3 x. |</>o(:x)| 2 = N. As usual in quantum 
optics, an Hamiltonian term involving a single quantum 
field operator leads to a coherent excitation of the field, 
in our case a collective excitation of the condensate. As 
expected, this term disappears in the case of a homo- 
geneous condensate considered in the rest of the paper 
thanks to the spatial orthogonality of the A(x) operator 
to the condensate wavefunction 0o( x )- 



d 3 x^(x) A(x) = 0. 



(20) 



In the general case, the density fluctuation pattern can 
be isolated even in the presence of a strong collective 
excitation simply by subtracting out the deterministic 
component of the density modulation. 



III. ADIABATIC LIMIT 



The mixing of the d k and dl k operators is negligible 
|Ck.+ 1 3> |Ck,- 1 — in the so-called adiabatic limit where 
the time-evolution a(t) tak es p lace on a very slow time- 
scale as compared to w k [45[ . As a consequence, the 
occupation rj k of the Bogoliubov modes is constant in 
time and equal to its initial value ?T- k h ' and the anomalous 
averages Ak remains zero. 

For a zero initial temperature T = 0, the adiabatic con- 
dition is equivalent to stating that the evolution is slow 
enough for the system to remain in its ground state: at all 
times, the density correlation function exactly coincides 
with the static T = one for the instantaneous value of 
a(t). On the other hand, for a non-zero initial tempera- 
ture T > even an adiabatic evolution is able to bring the 
system outside thermal equilibrium [46j | : while the pop- 
ulation n k of each Bogoliubov mode is conserved during 
the adiabatic evolution, the instantaneous energy ?kj k (t) 
of the mode has in fact a non-trivial time-dependence. 
As a result, the population n k of the different Bogoli- 
ubov modes at late times is no longer described by a 
simple thermal condition of the form (fT4)l . 



IV. SUDDEN JUMP 



General formulas 



Simple expressions for the expectation values appear- 
ing in (|18|) can be obtained in the limit of a sud- 
den variation of a(t) from a\ to a-i on a time scale 
at much faster than the frequency w k of all the rele- 
vant modes, i.e. fii t 2&t *C 1. In this limit of a sud- 
den quench [H, HE US EH, the evolution of the A k (£), 
A_ k (i) atomic operators during the sudden modulation 
of a is negligible and the simple picture of the dynam- 
ical Casimir effect introduced in [Tj can be straightfor- 
wardly applied: the d k , a_ k operators at t = 0^ right 
before and right after the jump are expanded in terms 
of the A k (i = 0), Al k (t = 0) atomic operators using ([5]) 
with the suitable u k , i> k Bogoliubov coefficients. An ex- 
plicit relation between the d k , d^ k operators at t = 0^ is 
straightforwardly obtained by eliminating the A k (t = 0), 
A_ k (i = 0) operators. The evolution of the a k , d_ k op- 
erators at later t > times reduces to the simple phase 
factor exp(=piw k t). 

Once all these steps arc combined together, one is fi- 
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nally led to the following compact relations [4l| : 
" 7 ?k + V k 



a k (t) 



■6k(0-) 



^M+_ k (cn) 



(21) 



? /k 



^k(o-) 



(o- 



where the ?7 k coefficients are defined as 
(w k /w k and the Bogoliubov frequencies lo 



(22) 
T before 



k / k / emu uiic uugmiuuuv iibqumubD u^ k 

and after the jump are evaluated using ^ with a = ai,2- 
As r]^ — > 1 for £k/Mi,2 3> 1, the transformations pTi 
do not mix the a k , al k operators for large values of k, 
which provides an a posteriori justification for the sudden 
jump condition o~t ^1,2 ^ 1- For the hydrodynamic modes 
with i? k *C Mi, 2, the mixing coefficient 7/ k tends instead 
to a finite limiting value (c2/ci) ±1 / 2 . 

The Bogoliubov mode occupation n k after the sudden 
change of scattering length is given by the formula 



+ ( w k ) 2 th,0 



fa 



2uj+ w k 



4w k w k 



while the anomalous average has the form 



2n 



th.O 



1 



-2zw, + t 



(23) 



(24) 



Note that the Bose distribution n^'° is here to be eval- 
uated according to (|14[) using the initial value w k of the 
Bogoliubov mode frequency. In agreement with (jl~3)) . the 
initial value of the anomalous average has been taken to 
be zero. The terms in (|2"3"|) and proportional to the 
initial population n, k h '° account for the amplification of 
initial thermal excitations by the sudden jump in a, while 
the other terms describe the contribution due to the dy- 
namical Casimir emission of Bogoliubov phonons out of 
the initial vacuum state. 

It is interesting to evaluate (f2"3")) and in the limit of 
a small change of a(t), i.e. 02 = 0,1 + Aa with Aa <C 01. 



In this limit, one has: 



"k 



th.O 



A/1 2 



4(£ k + 2fj,y 



A k (t) 



A /Li 



2{E k + 2/i) 



2// 



th,0 



2n 



1 



th,0 



-2iw, + t 



e 



(25) 
(26) 



Here, the chemical potential variation is defined as A/i = 
Anti 2 (a2 — ai)no/m. While the effect on the anoma- 
lous average is linear in A//, the population change is 
quadratic and therefore much weaker. The difference 
is even more dramatic at T > 0: while the population 
change is a (small) correction proportional to A/1 2 on 
top of the (large) initial thermal distribution, the anoma- 
lous average fully originates from the dynamical Casimir 
emission and is amplified by the initial thermal popula- 
tion. 



B. Physical discussion 

This physics is illustrated in the plots of the Bogoli- 
ubov occupation and the modulus of the anomalous av- 
erage shown in Fig[TJd,e) for the 0,2/ a\ = 0.25 case. For 
these plot, the exact formulas ([23)) and (|24|) have been 
used. 

For a vanishing initial temperature T = (black line), 
both n k and „4k show a smooth peak centered at k = 
and a power-law tail that extends far on the high en- 
ergy modes. As predicted by the analytical approximate 
formulas (|2"5")) and (|2l)|) . the anomalous average *4 k is in 
modulus much larger than the occupation 77, k . Note that 
the smooth peak would be replaced by a 1/fc divergence 
if interactions in either the initial or final states were 
vanishing. 

For a finite initial temperature T > 0, the effect of the 
sudden variation of a on n k is a very weak correction on 
top of the initial thermal distribution n k h '° '. On the other 
hand, the T = contribution to the anomalous average 
due to the dynamical Casimir effect is strongly amplified 
by the initial thermal populaton. 
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FIG. 1: Panels (a-c): Time evolution of the density corre- 
lation function after a sudden jump of the scattering length 
from ai to 02 = ai/4. Different (a-c) panels refer to dif- 
ferent evolution times after the jump, fi\t = 0, 5, 20. Panel 

(d) : Bogoliubov mode occupation n k after the jump. Panel 

(e) : anomalous average *4 k | after the jump. In all panels, 
thick black lines correspond to an initial T = 0; thin red lines 
correspond to a finite initial temperature T//11 = 0.5. The 
dashed red line in panel (d) indicates the prediction of the 
adiabatic model, i.e. the initial population of the Bogoliubov 
mode n k h '°. 



The density correlation function is immediately ob- 
tained inserting (|23ti24|) into the general expression (| 18|) . 
A few snapshots after different evolution times arc shown 
in FigQJa-c) for the simplest case of a one-dimensional 
system. Exception made for some geometrical factors, 
the physics is however identical when two- or three- 
dimensional systems are considered. 

For t — + , straightforward algebraic manipulations 
confirm that the sudden change of a does not have an 
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immediate effect on correlation function, g^(x— x!;t — 
+ ) = 5 (2) ( x - x';< = 0"): The jump is in fact too 
rapid for the microscopic state of the atoms to respond. 
However, as this state is no longer an eigenstate of the 
system Hamiltonian with a(t > 0) = a 2 , a non-trivial 
evolution is observed on g( 2 ' at later times t > 0. 

Before the sudden change of a, the density correlation 
function is characterized by a dip around x = x' due to 
the effect of atom-atom repulsive interactions. At finite 
temperature, this dip is less pronounced than at T = 0, 
and starts being accompanied by a Hanbury-Brown and 
Twiss bump due to the thermal fluctuations 

After the sudden change, the static central structure 
around x = x' is somehow amplified by the change in 
population and, more importantly, by the increase in 
the Bogoliubov + coefficient as a consequence of 
the reduced value of a. At the same time, a system of 
moving fringes originates from x = x' and propagates 
in the outwards direction as a consequence of the time- 
dependent anomalous average Ak(t). At each time t, the 
fringe pattern is concentrated in the \x — x'\ > c 2 t region 
and shows a significant chirping in space, the external 
part of the pattern having a shorter wavelength than the 
inner part. As time goes on, the fringe pattern gets pro- 
gressively stretched in space. 

This peculiar fringe pattern has a very simple physical 
interpretation [25[. When the sudden change of a occurs, 
pairs of entangled phonons are created at all spatial posi- 
tions with opposite momenta ±fc. As time goes on, these 
pairs propagate in opposite direction at a fc-dependent 
group velocity v 2 ,k- As the group velocity v 2 ^ of Bogoli- 
ubov modes is a growing function of k which starts from 
V2,k=Q = C2, the correlated pairs will be separated at a 
time t by a distance equal or larger than 2c 2 t [47| More 
precisely, the modes that most contribute to the fringe 
pattern for a given \x — x'\ are the ones with a wavevector 
±/c such that v 2 j — \x — x'\/2t, which are responsible for 
fringes of wavelength 1/fc. In this semiclassical picture, 
the observed chirping is then a simple conseguence of the 
fact that higher k modes propagate at a faster group ve- 
locity. 

While the strongly chirped external region of the fringe 
pattern remains almost unaffected by a finite initial tem- 
perature, the long wavelength fringes at low \x — x'\ are 
substantially reinforced. This confirms our intuitive un- 
derstanding of the fringe pattern: according to (|2~i)) ther- 
mal enhancement is in fact concentrated into the low- 
k modes which are responsible for the long wavelength 
fringes. 



V. SLOW RAMP 

The previous Section was devoted to the case of a sud- 
den change of a for which analytical expressions were 
available for the Bogoliubov mode amplitudes in the final 
state. The more general case of an arbitrary dependence 
a(t) requires a solution of the pair of ordinary differential 



equations (|9l fT0|) . In the present section we discuss the 
result of a numerical solution of these equations for the 
case of a smooth temporal dependence of the Erf form: 

ai + ft9 a 2 — ai ft — to \ 
( t ) = ^ + JL_LErf(— °J. (27) 

where the change of scattering length from a\ to a 2 takes 
place on a time scale at- 




FIG. 2: Panels (a-c): Time evolution of the density cor- 
relation function after a slow ramp of the scattering length 
from ai to a,2 = ai/4. The ramp follows an Erf shape with 
a t = 5/ fii centered at t — 5a t . Different (a-c) panels refer 
to different evolution times after the jump, fiit — 0, 50, 70. 
Panel (d): Bogoliubov mode occupation after the jump. 
A magnified view of the peak is given in the inset. Panel (e): 
anomalous average |.4k| after the jump. In all panels, thick 
black lines correspond to an initial T = 0; thin red lines corre- 
spond to a finite initial temperature T/fil = 0.5. The dashed 
red line in (d) is the initial population n^'° . 

As one can see in FigJ2£d,e), the main effect of a finite 
a t is to introduce a ultraviolet cut-off to the Bogoliubov 
modes that are effectively excited during the modulation 
of a. All Bogoliubov modes with u>ka t 3> 1 experience 
in fact the modulation as adiabatic: as a consequence, 
the corresponding anomalous average Ak remains fully 
negligible and the population remains very close to 
its value before the ramp [dashed line in Fig[2fd)]. 

The density correlation function at different times af- 
ter the slow modulation is shown in FigJ21(a-c) for the 
simplest case of a slow modulation rate as compared to 
the chemical potential, fj-i^t 2> 1- In this case, only the 
low-fc, hydrodynamic modes result appreciably excited 
and the chirped fringes in the large \x — x'\ 3> 2c 2 t region 
disappear from the fringe patterns shown in Fig[2Ja-c). 
These are then characterized by a single negative peak 
that rigidly propagates at a speed 2c 2 with almost no 
dispersion. A different point of view on this same phe- 
nomenology was recently presented in [25[. As one can 
see by comparing the thick black lines to the thin red ones 
in Fig[^a-c), in this case the effect of a finite initial tem- 
perature reduces to an amplification of the propagating 
peak. 
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VI. HYDRODYNAMIC LIMIT 

In the limit of a slow ^iat/h 1 and weak |A/x| <C n 
jump, an analytical approximation can be obtained for 
the height and shape of the moving peak. The idea is 
to use the sudden jump result (|26|) and then take into 
account the slow variation of a(t) by means of a cut- 
off in the momentum integration of (|18j) : while the low- 
frequency modes Wk <C l/<r t experience the modulation 
as sudden, the high-frequency ones lu^ 1/ct experi- 
ence it as adiabatic. The assumed condition fj,at/h 3> 1 
implies that the momentum cut-off is at a wavevector 
fc^x ~ 1/coj <C l/£ w cll within the hydrodynamical 
regime for which ~ c |k|. 

Including this cut-off as an additional exponential fac- 
tor exp(— fc/fc^ax) m the integral of the zero-point con- 
tribution to (US]), we immediately get to the following 
expression for the moving peaks: 



Depending on whether k B T ^ , the height of the 



thermal peaks is either: 



x ( 2 ) 
S 9th 



1 



An k B T 



peak,high-T 4 (n£) (flUt/H) fi 2 



X ( 2 ) 
^th 



1 Afi (kBT) 12 



poak,low-T 4(n£) 



(31) 



(32) 



Even though the parameters used in FigJ5] are on the 
edge of the validity domain of the hydrodynamic approx- 
imation, the analytical formulas discussed in the present 
section turns out to be in reasonable quantitative agree- 
ment with the numerical results. 



VII. COMPARISON WITH BLACK HOLE 
CALCULATIONS 



!>9t=o( x > x ') 



1\-{x 



' - 2ct) 2 



(x 



2ct) 2 



[£ 2 + {x - x> - 2ct) 2 } 2 \l\ + (x - x> + 2ct) 



212 



The peak value is at 



x ( 2 ) 
°9t=o 



i An 

peak 47T (n£) {\MJt / K) 2 M 



(28) 



(29) 



and the width is determined by the cut-off as it = 
1/^max = ca t '- Physically, this value of the peak width 
can be understood as a consequence of the uncertainty 
a t in the emission time, which reflects into a broadening 
It = cat of the correlation signal. 

At a finite temperature T > 0, one has to include 
the further contribution due to the amplified thermal 
fluctuations. Depending on whether the temperature 
k B T is higher or lower than the effective cut-off energy 
-^max = ^Lax = ft/of imposed by the slow ramp, the 
cut-off on the thermal contribution to (fl"8)) has to be im- 
posed at A:JJ ax = min(/c^ lax , k B T/Hc). Including again 
this cut-off as an exponential factor exp(— fc/fc^ ax ), one 
gets the following expression for the moving peaks: 



<*ffth \ x i x ') 



Afi k B T i 



th 



2/i 2 n 



1 



(x - x> - 2ct) 2 + e 2 h 



(x 



2ct) 2 + 1 



(30) 



The width of the thermal peaks is set by the cut-off £th = 
1/^max- At l° w temperature k B T < E^ laK , the width 
is enlarged to £th = frc/k B T as a consequence of the 
finite correlation length of thermal density fluctuations 
in the initial state. At high temperature k B T > E^^, 
the width is again dominated by the ramp time effect as 
in the T = case. 




40 60 

\x-x']/L 



FIG. 3: Comparison of the Bogoliubov prediction (|18[1 for the 
density correlation function (thick black line) with the result 
of the numerical simulations of [22| (thin red dashed line). 
The numerical lines are cuts of G^ 2 \x,x') along a x + x = 
2io line with xo/£i,2 2> 1 well inside the acoustic black hole 
region; all the field of view is within the acoustic black hole 
region. The scattering length is brought from ai to a,2 = ai /4 
with a Arctan temporal dependence on a time scale a t /-ii = 
0.5. The correlation functions are evaluated a time flit = 50 
after the change of a. The upper (a) panel is for T = 0. The 
lower (b) panel is for a finite ksT j \i\ — 0.1. 



In the previous sections we have investigated in detail 
the density correlations that appear as a consequence of 
a modulation of the atomic scattering length a(t). One 
of the motivations of the present work was to fully un- 
derstand some unexpected transient features that were 
observed in the numerical experiment of [221 ] . namely a 
system of moving fringes that appear inside the acous- 
tic black hole as soon as the horizon is created and then 
rapidly leave the field of view. As the acoustic horizon 
was created by suddenly ramping down the atomic scat- 
tering length in a full half space and a quantitatively 
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identical system of fringes was observed in a spatially 
homogeneous system, an interpretation was put forward 
in terms of dynamic Casimir effect. In this Section, we 
confirm this interpretation by performing a quantitative 
comparison of the numerical results of [22| to the predic- 
tions of the Bogoliubov model that we have discussed in 
the previous Sections. To this purpose, the same Arctan- 
shaped ramp of a(t) that was used in the numerical cal- 
culations has to be implemented in the Bogoliubov calcu- 
lation: the results of the comparison are shown in FigJ3] 
The agreement between the two calculations is remark- 
able, which firmly confirms our initial interpretation. 



VIII. QUASI-PERIODIC MODULATION 




k B T/n 



FIG. 4: Panels (a,b): Time evolution of the density correla- 
tion function after an oscillating modulation of the scattering 
length of amplitude 8a/ a — 0.1, carrier frequency Tiwo/a 4 = 1, 
and Gaussian envelope of duration uqT = 10. The two panels 
refer to different evolution times after the peak of the modu- 
lation, fj,t = 150 (a) and 200 (b). Panel (c): Bogoliubov mode 
occupation rik after the modulation. Panel (d): anomalous 
average |„4k| after the modulation. In all (a-d) panels, black 
lines are for a zero initial temperature T — 0; redlines are for 
an initial temperature ksT/fi = 0.5. Panel (e): Peak value of 
the fringe amplitude as a function of the initial temperature 
for a given quasi-periodic excitation sequence. Black circles: 
numerical integration of (|9I10[) . Dashed red line: fit of the 
points with a thermal law (1 + 2nJ t h '°). 



A narrow window of k modes can be specifically ad- 
dressed by using a periodic modulation of the scattering 
length in time: according to the form (JSJ) of the system 
Hamiltonian, a weak perturbation at frequency luq is in 
fact able to effectively excite those pairs of Bogoliubov 
modes that satisfy the resonance condition 



LU = W k + W_ k - 



(33) 



This physics is illustrated in FigJH where we show the 
result of a numerical integration of ([51 ll0p under a sinu- 
soidal modulation of a(t) with a Gaussian temporal en- 



velope: the mixing of the <2k and a^L k operators is limited 
to a small range of k vectors and results in very peaked 
shapes of the Bogoliubov mode occupation and of the 
anomalous average^ [panels (c,d)]. 

The pair of weaker peaks that appears in the anoma- 
lous average at larger values of k is due to second- 
order processes in the modulation amplitude. This inter- 
pretation is confirmed by the scaling of the peak ampli- 
tude as 5a 2 and by the position of the peak that satisfies 
the second-order resonance condition Uk + u;_k = 2ujq. 
For larger modulation amplitudes, higher order peaks 
would appear at k values satisfying higher-order reso- 
nance conditions of the form lu^ + w_k = N u>q, N being 
a generic positive integer number. 

Correspondingly to the dominant pair of peaks, the 
density correlation function [panels (a,b)] shows a peri- 
odic fringe pattern of wavelength 2ir/k that travels away 
from \x — x'\ = at the group velocity 2v2,k- The en- 
velope of the fringe pattern follows the envelope of the 
oscillating a(t) modulation. 

At a finite temperature (thin red lines) , the resonance 
peaks in the Bogoliubov mode occupation are almost 
completely hidden by the thermal component, but re- 
main perfectly visible and even reinforced in the anoma- 
lous average A^. Correspondingly, the moving fringe 
pattern experiences an overall amplification without any 
significant distortion of the oscillating shape nor of its 
envelope. 

The simple relations (fTS)) and ([To]) that relate the den- 
sity correlation at the end of the modulation sequence 
to the initial thermal occupation of the mode can be ex- 
ploited as a simple way to precisely measure the tem- 
perature of the system in an almost non-destructive way. 
This proposal extends an original suggestion of [24[ to 
measure the temperature of a Bose-Einstein condensate 
using a parametric modulation of some parameter: as il- 
lustrated in FigfJJc), looking at the density correlation 
rather than at the Bogoliubov mode occupation has the 
significant advantage that the interesting signal is not 
hidden by a broad thermal pedestal. In contrast to the 
case of a single jump discussed [28j |. a periodic modula- 
tion is able to concentrate the interesting signal into a 
single Bogoliubov mode and, more remarkably, to make 
it significantly stronger without distorting it. 

As long as the applied modulation is weak enough for 
nonlinear and saturation effect beyond Bogoliubov the- 
ory to be negligible, the observed signal is in fact propor- 
tional to 2nj c h '° + 1. To clarify this statement, the peak 
value of the fringe amplitude is plotted in FigHJe) as a 
function of the initial temperature. The points are the 
result of a numerical integration of ([ WTII)) . the dashed line 
is a fit using the known thermal dependence: provided a 
suitably low-energy mode is used, the peak value of the 
fringe amplitude is proportional to the system tempera- 
ture. 
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IX. REINFORCING THE DENSITY 
CORRELATION 

A critical issue in view of an experimental verifica- 
tion of the conclusions of the present paper as well as 
of the predictions of Q is the actual value of the den- 
sity correlation signal that one is to detect: given its 
relatively small value, methods to reinforce it can be of 
crucial importance. In the present section we apply to the 
dynamical Casimir effect a diagnostic trick that was re- 
cently used to characterize phase fluctuations of a quasi- 
condensate in a strongly one-dimensional geometry e.g. 
in [43l | and that was recently put forward in the context 
of observing the analog Hawking radiation [35[ . The effi- 
ciency of this tool to measure the microscopic properties 
of low-dimensional many-body systems was recently dis- 
cussed in [28[. A related idea was presented in [l7| with 
the purpose of amplifying the signal of analog cosmolog- 
ical particle production. 
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I x-x' I / 4 

FIG. 5: Density correlation function after a slow ramp of the 
scattering length from a± to 0,2 = ai/4 as in Fig(5] a faster 
switch-off of the scattering length to aj = within a' t <C 
a t , and a final time interval t{ ICC of ballistic, non-interacting 
evolution. Different panels (a-c) correspond to different free 
evolution times fii tf lac = (a) 5 (b), and 10 (c). The switch- 
off of the scattering length to a/ = is performed at jj.it = 70 
within a time fJ,i<r' t = 1 (black lines). For comparison, the 
case of a sudden switch-off (a' t = 0) is shown as a red dotted 
line in (b). 



In agreement to the Goldstone theorem, long wave- 
length phonons have a mostly phase-like character and a 
very weak component of density fluctuation (37j : a quick 
switch-off of the the interactions shortly before measur- 
ing the density correlations can then be used in order to 
reinforce the signal by converting phase fluctuations into 
density fluctuations. The efficiency of this trick is illus- 
trated in FigO the density correlation signal is plotted 
for different values of the ballistic expansion time ts n be- 
tween the switch-off of a and the actual measurement. 
During this time, the original signal gets amplified by a 
significative factor. 



Note that in order to avoid a substantial emission of 
high-fc particles and the consequent appearance of fast 
oscillations in the density correlation pattern, the switch- 
off time <j' t can not be chosen too short. This point is 
illustrated in Figj5^b) where the signal obtained with a 
sudden final jump is shown for comparison as a red dotted 
line: the importance of a careful choice of a' t is apparent. 

An analytical understanding of the physical origin of 
the different features that appear in the density correla- 
tion function after the second jump is the subject of the 
next subsection. 



A. Hydrodynamic model 



Analytical expressions can be obtained in the case in 
which the second jump brings the scattering length to a 
finite final value a/ and both jumps are performed on a 
time-scale at , o' t long as compared to the chemical poten- 
tial, C7t,u' t 3> Under these assumptions, the ana- 
lytical technique introduced in Sec lVll can be generalized 
to the case of a two-jump modulation sequence. 
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FIG. 6: Density correlation function after a two-jump se- 
quence. Jump times at, a' t = 8/fii. Delay time tdei = 100//xi. 
Upper panel (a): 02 = ai/4, a/ = ai, observation time 
t — ffroc + tdd = 270//zi. Lower panel (b): 02 = ai/4, 
a/ = ai/8, observation time t = tf Ica + tdei = 320//ii. The 
blue labels indicate the terms in the two-sudden-jumps ana- 
lytical model that correspond to each feature. Red lines indi- 
cate the density correlation function in the absence of second 
jump, af = ai- 



An explicit forms for the Bogoliubov coefficients a time 
£fi n after the end of the two-jump modulation sequence 
can be obtained by repeatedly applying two transforma- 



10 




suit: 



FIG. 7: Scheme of the spatial position of the different features 
predicted by (|38p for a configuration ai > a2 > a/ inspired 
to FigEIb). 



tions of the form (I21II22I) 



C, 



k,+ 



e k 



r (2) r (l) _-i4 2 ' W 
(2) ^(1) O i4 2) t dol 



k.+ ' 



(34) 



(35) 



Here, td e \ is the time interval between the jumps and ta n 
is the time interval between the second jump and the 
actual measurement, uj^ ' are the Bogoliubov disper- 
sions before the first jump, after the first jump, and after 
second jump, respectively. The single jump Bogoliubov 
coefficients C^'^ are defined according to (|21|) and 
as: 



(i) 

k,± 



°k,± - 



i 



4 2) 
4 1} 




4 X) 
4 2) 


4 /} 
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4 2) 


4 2) 


4 /} 



(36) 
(37) 



The corresponding density correlation function is then 
obtained by inserting these formulas into the general for- 
mula (|18p and imposing a suitable cut-off to the integrals 
at fcJnax = l/£t- Limiting ourselves to the simplest T = 
case, some straightforward algebra leads to the final re- 



5g^{X = x-x') 



{A x F it [X - 2c 2 t del \- 



2Trmncf 

+ A Ft, [X] + A- X F u [X + 2c 2 t dcl } + 
+ Bi F lt [X - 2cztft.ee - 2c 2 t d ei] + 

+ B Fg t [X — 2c/tf ree ] 
+ S_! F tt [X - 2 C/ iftee + 2c 2 t d el] + 

+(X*->-X)}. (38) 

The amplitudes have the following expressions in terms 
of the Bogoliubov operators of the two jumps: 

(39) 
(40) 
(41) 

(42) 
(43) 
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and the function F/(x) is defined as 
Fi{x) 



(44) 



As a consequence of the interference between the dif- 
ferent terms of (|34H35[) , a number of peak/dips appear in 
the final result ([38]) and have a peculiar evolution as a 
function of tf ree . An illustration of this physics is shown 
in FigO even though significantly distorted by effects 
beyond hydrodynamics, all the features are clearly recog- 
nizable. Labels refer to the corresponding amplitudes de- 
fined in eqs. (|39H4"3")) and schematically illustrated in Fig|7] 

The Aq feature provides a slight modification of the 
many-body dip. The standard dynamical Casimir effect 
by the second jump is responsible for the feature Bq that 
emerges from the many-body dip at x = x' at travels at 
a speed 2c/. In agreement with (|29[) . its sign depends on 
the sign of the second jump Aaf = af — a 2 - 

The dynamical Casimir feature that was visible at 
2c2idei before the second jump splits into three features 
B—i, A±i and B\ that travel away at speeds respectively 
equal to — 2c/, 0, 2c/. In the absence of second jump (i.e. 
for af — ct2, dashed red line in FigO, only the B\ sur- 
vives with a finite amplitude, though at a slightly shifted 
position as a consequence of the unchanged sound veloc- 
ity. All other -B_i, A±\ features instead vanish as a con- 

(2) 

sequence of the e^ _ = condition. Before the jump, the 

height of the B\ feature is proportional to Cg_/c2. 
At the jump, it gets multiplied by a factor approximately 
equal to: 



£2 



(45) 
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As expected, this factor is larger than 1 as soon as the 
second jump corresponds to a decrease in the scattering 
length af < a,2- In particular, it becomes very large when 
a.f is brought to a very small value af <C a 2 ■ 

Even though these analytical considerations are lim- 
ited to the hydrodynamic regimes, they provide an use- 
ful qualitative guidelines to interpret the full numerical 
results shown in Fig [5] and [|5] 

X. CONCLUSIONS 

In this paper, we have presented a general theory of 
the density fluctuations that appear in an atomic Bosc- 
Einstein condensate as a consequence of a temporal mod- 
ulation of the atomic scattering length. Different regimes 
have been identified as a function of the time scale and 



the temporal shape of the modulation. A physical pic- 
ture in terms of the dynamical Casimir emission of pairs 
of entangled phonons has provided an intuitive explana- 
tion of the results. Simple analytical formulas have been 
obtained in the most remarkable limiting cases. Excel- 
lent agreement is found with the quantum Monte Carlo 
calculation of Ref. [23], which a posteriori confirms the 
physical interpretation of the numerical data. The effi- 
ciency of a recently proposed strategy [35[ to reinforce 
the experimental signal is discussed and quantitatively 
validated. Possible applications to the thermometry of 
ultracold atomic gases are pointed out. 

We are grateful to C. Tozzo, F. Dalfovo, E. Cornell, 
and P. Calabrese for stimulating exchanges and discus- 
sions. A long-lasting collaboration with C. Ciuti and S. 
De Liberato on the Dynamical Casimir Effect is warmly 
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